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Abstract 

We consider the implications of solar matter density random noise 
, upon resonant neutrino conversion. The evolution equation describing 

I MSW-like conversion is derived in the framework of the Schrodinger 

■ approach. We study quantitatively their effect upon both large and 

, small mixing angle MSW solutions to the solar neutrino problem. This 

\ is carried out both for the active-active Ve i^fi,T as well as active- 

' sterile fg fg conversion channels. We find that the small mixing 

. MSW solution is much more stable (especially in Am^) than the large 

\ mixing solution. The possible existence of solar matter density noise 

at the few percent level could be tested at future solar neutrino ex- 



Ph. periments, especially Borexino. 
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1. Introduction 



The long-standing deficit of solar neutrinos (the Solar Neutrino Problem 
(SNP)) has now been observed by all four operating experiments . 
The main essence of the SNP is the strong deficit of the beryllium neutrinos 
1^. On the other hand, the high energy boron neutrinos are moderately 
suppressed, while the low energy ones are almost undepleted. This strongly 
suggests that any astrophysical solution fails [^,0 in reconciling the experi- 
mental data with the Standard Solar Model (SSM) predictions [p|,p|,p!0|. 

It is possible to ascribe the solar neutrino deficit to the existence of two 
types of neutrino conversion mechanisms, both of which can deplete neutrinos 
of different energies differently, as required by the experimental data. The 
long wavelength vacuum oscillations provide a good fit to the most recent 
results for Am^ ~ 10~^°eV^ and large neutrino mixing sin^ 26* ~ 1 [|lT],0,|l2 



The other scenario is the resonant neutrino conversion due to interactions 
with constituents of the solar material (the Mikheyev-Smirnov-Wolfenstein 
(MSW) effect) |]14[. This provides an extremely good data fit in the small 
mixing region with Am^ ~ 10~^eV^ and sin^ 20 ~ 10~^ 10~^ [|15|,0,|12 



Both of these solutions have been studied against possible changes of the 
SSM input parameters For example, the study of the MSW effect 

has revealed its stability, especially in the Am? parameter. 

In this paper we investigate the stability of the MSW solution with respect 
to the possible presence of random perturbations in the solar matter density, 
so far not included in the standard MSW picture. 

In Ref. [0 the effect of periodic matter density perturbations added to 



an average density po? i-e. 

p(r) = po[l + /isin(7r)] (1.1) 

upon resonant neutrino conversion was investigated. The major effects show 
up when the fixed frequency (7) of the perturbation is close to the neutrino 
oscillation eigen- frequency, and for rather large amplitude values (/;, ~ 0.1 — 



0.2) , giving rise to the parametric effects |T^ . Such effects can either enhance 
or suppress neutrino conversion in the Sun. There are also a number of papers 
which address similar effects by different approaches |[T^JT^. 



Direct observations of solar surface motions, resulting from the superpo- 
sition of several modes, may indicate a rich spectrum of frequencies. This 
would suggest the need to consider the effect of random or "white" noise mat- 
ter density perturbations ^(r), characterised by an arbitrary wave number 

^(r) = J dk^{k)smkr , (1.2) 

rather than a periodic or regular perturbation. In such the spatial 

correlation function for a uniform medium 

(e(ri)^(r2)) = {e)r.-r„ (1.3) 

obeys {C{k)^{k')) = {^^)kS{k + k') as the averaging rule for the Fourier com- 
ponents, where the wave number k is not fixed. The effect of solar density 
as well as solar magnetic field fluctuations upon neutrino spin-flavour con- 
versions has also been considered in Ref. |T^, using somewhat different 
methods. 



In this paper, after some discussion (Sec. 2) about the nature of the 
matter density fluctuations, we derive the most general neutrino evolution 
equation in random matter, starting from the standard Schrodinger equation 
(Sec. 3). This discussion is closer to the particle physics intuition than that 
of Ref. |[19|. Moreover, we consider both the active-active z/g — u^ j. as 
well as the active-sterile z/g — >■ z/^ neutrino conversion channels (here Ug is a 
neutrino state with no standard model interaction). The latter is motivated 
by the fact that the existence of a sterile neutrino seems to be the only way 
to simultaneously account for the solar and atmospheric neutrino deficits in 



the presence of neutrino dark matter 



After an analytical study of the neutrino conversion equations we have 
investigated the impact of matter density noise upon the MSW scenario in 
the context of the SNP (Sec. 4). Typically, we find that the presence of mat- 
ter fluctuations weakens the MSW mechanism, thus reducing the resonant 
conversion probabilities |T9|. We have carried out a fit of the latest solar 



neutrino data for different values of the noise level, minimising the the 
(Am^, sin^ 26) plane. As in the noiseless MSW case, we find that the small 
mixing MSW solution provides a better fit to the data than the large mixing 
one, both for the case of active, as well as sterile neutrino conversions. We 
present several plots with the results of our fits in which the effect of the noise 
is studied in the idealised approximation where all neutrinos are produced 
at the solar centre. We conclude that the most relevant parameter region 
corresponding to adiabatic conversion of ''Be neutrinos is relatively stable 
with respect to such density fluctuations, whereas there is a larger effect of 
the noise for the large mixing MSW solution. We show how the possible 
existence of solar matter density noise could be tested in the next generation 
of solar neutrino experiments, especially Borexino. Finally, we comment on 
how possible solar model uncertainties could affect our results. 



2. Matter Density Noise in the Sun 

Let us briefly discuss the expected size of fluctuations in the Sun and their 
correlation lengths. For the sake of discussion, we can approximate (except 
in the very inner core) the average solar matter density, as given by the 
Standard Solar Model (SSM) [|,|,|Tg, by: 

p(r)~poexp( ^) (2.1) 

/to 

where Rq ~ O.lRs {Rg is the solar radius), tq — 0, and po ~ 250g/cm^. The 
SSM in itself cannot account for the existence of density perturbations, since 
it is based on hydrostatic evolution equations. 

One may however speculate upon possible mechanisms that could induce 
such density inhomogeneities in the Sun. Unfortunately it is quite difficult 
to give reliable estimates of the density perturbations in deep layers of the 
Sun, since this would require the detection of (^-modes, not yet possible [PT|] . 
Indeed, these modes can exist only in deep layers beneath the convective 
zone and thus they can reach the surface only after an exponential damping 
through the convective zone 0. Note also that it is extremely difficult to 
identify the ^f-modes in helioseismology observations, due to their tendency 



to be accumulated in the lower frequency part of the Fourier spectrum. We 
may, however, give a simple estimate of the level of density perturbations 
6p in the solar interior by combining the continuity equation up to the first 
order in 5p and the velocity perturbation 5v 

^ + V ■ 5vp = 0, (2.2) 

with the p-mode observations of the IRIS network at Tenerife 0. These 
show that in the lower frequency part of the Fourier spectrum, the p-mode 
spectrum resembles that of noise, namely ~ 1//. For instance, using the 
measured power 5P = lO^m^s^^ corresponding to the frequency / = 10~^ 
s~^ (see Fig. 26 of Ref. 0) we may estimate from Eq. ( p.2|) the perturbation 
level ^ defined as 



« ^ ^ . (2.3) 

P P 



We obtain 



^ ~ PLl fLl ^ lO-^s-i X (106m)2 ^ ^ - ' 

where we have taken as typical size of the spatial inhomogeneity the value 
Lq ~ 10^ km, the so-called "granule" -size. Thus we see that values ^ ~ 0.3% 
in the solar surface can not be excluded. In contrast, inside the solar core 
the estimate of the parameter Lq becomes very rough. In fact one expects 
that, due to buoyancy, the (7-mode amplitudes beneath the convective zone 
can be larger than at the surface and, correspondingly, the inhomogeneity 
size Lq smaller than at the edge of the Sun (see Fig. 13 (a) of Ref. 0). As 
a result, for a fixed perturbation amplitude 5f ^ = const we can extrapolate 
the continuous power spectrum to low frequencies leading to a large density 
inhomogeneities since ^ ~ /"^Lq ^. 

There is another way to estimate the level of density fluctuations using 
the density profile of Eq. (|2.1|) . Indeed, in the hydro-dynamical approxima- 
tion, density perturbations can be induced by corresponding temperature T 
fluctuations due to convection of matter between layers with different local 
temperatures. For example, if we express the macroscopic matter density p(r) 
through the Boltzmann distribution with the gravitational potential energy 
U = mpg{r){r — tq), where g{r) = GM{r)/r'^, G being Newton's constant, nip 
the nucleon mass and M(r) the mass contained in a sphere of radius r. The 
change T — >• T + (5T leads to 

p(r ; T + 6T) = po exp [- ^^^f^^^] = p(r ; T) [1 + 

where ^ = C,{r; T; 6T). From this we have 

^ _ Sp _ nipgjr - rp) 6T _ {r - rp) 5T 

^ p T T i?p T ' ^ ■ ^ 

where we have compared the relevant exponent with that in eq. (|2.1j). One 



can argue that J {6T'^) /T < 0.05 is not in conflict neither with the SSM 



nor with present helioseismology observations H,22]. Thus from Eq 



taking [r — tq) / Rq ~ 1, we obtain a correspondingly comparable level of 
density fluctuations. Thus in what follows we assume the existence of such 
few percent level matter density fluctuations, up to 8% Q. 

Now we generalise the above discussion to the case in which the perturba- 



tion 5p is of random nature. Following ref. |T9[ we assume that the random 
field 5p is a 5-correlated Gaussian distribution. For small inhomogeneities, 
the autocorrelation function (,^^) can be taken as 

(5p(ri)5p(r2)) = 2p'{e)L,5{n - r^) (2.6) 
whose correlation length Lq obeys the following relation: 

Ifree Lq <^ Xm (2.7) 

where /free = (cr^o)~^ is the mean free path of the electrons in the Sun. This 
lower bound is dictated by the hydro-dynamical approximation used later. 
For Coulomb interactions, the cross-section a is determined by the classical 
radius of electron roe = e'^/nieC^ ~ 2 x 10~^^cm, resulting in /^ee ~ 10 cm for 
a solar mean density uq ~ lO^^cm"^ and cr ~ 10~^^cm^. On the other hand, 
the upper bound expresses the fact that the scale of fluctuations should be 
much smaller than the characteristic neutrino matter oscillation length. Am, 
as indeed the 5-correlation distribution in Eq. (p.6|) requires. 



3. Neutrino Conversion in Noisy Matter 

Let us consider a system of two neutrinos and Ux. In the case of active- 
active neutrino conversion x = fi or t, while for the case of active-sterile 
neutrino conversions = Vsi with Vg being the sterile state. 

Neutrino propagation in the solar medium is affected by the coherent 
neutrino scattering off matter constituents which can be described in terms 
of the matter potential V . In the rest frame of the unpolarised matter, the 
potential is given, in the Standard Model, by 

V = ^^pY (3.1) 

nip 

where Gj? is the Fermi constant, p is the matter density and F is a num- 
ber which depends on the neutrino type and on the chemical content of the 
medium. More precisely, y = Fg — |y„ for the Ve state, Y = —\Yn for z/^ and 
Vr and y = for the state, where Fe,n denotes the electron and neutron 
number per nucleon. The matter potential modifies the energy dispersion 
relations for neutrino states, leading to the phenomenon of resonant conver- 
sion (the MSW effect 0]). Let us note that in this respect the potential V ^ 
i.e. the function p, previously described, represents an average macroscopic 
quantity. 



"'^Note that for an ideal plasma, like that in the Sun, the equilibrium plasma fluctuations 
are neeuE ible 11, i ^<5nl>^/^ /ue - Kr^)"^ < 1, where Ue is the electron density, 
Sue the corresponding fluctuation and — {T / 'iire^neY^'^ the Debye radius. Since the 
number of particles inside the Debye radius is very large, Nd = ner\) ^ 1, with ro ^ f 0"^ 
cm Ifree thcsc fluctuations are irrelevant for our present discussion. 



Now we re-derive the evolution equation for the neutrino in the presence 
of matter density random perturbations, which we regard as superimposed 
over the main average matter density profile. It is clear from Eq. ( p.6| ) that 
the random component of the potential can be written as V^(t)^. 

The evolution for the v^ — Vy {y = x or y = s) system is governed by 



d f f Hp H, 



'dAvy)-[Hly Hy)[uy)' ^^"^^ 

where the entries of the Hamiltonian matrix are given by f\ 

H, = 2[A,y{t) + A,y{t)], Hy=0, 

1 Am^ ~ 1 

Aeyit) = -[V^yit) -—^cos2ei A,y(t) = -Veym (3.3) 

Here 6 is the neutrino mixing angle in vacuum, Am^ is the neutrino squared 
mass difference, and the matter potential for the active-active neutrino conversion 
{y = x) reads 

K.(t) = ^^p(t)(i-v;) (3.4) 



rup 



or alternatively in case of 



VUt) = ^^pm-lYr,) (3.5) 



rUp 



(the neutral matter relation Yf, = 1 — Y^ has been used). 

The above system can be rewritten in terms of the following equations: 

P{t) = 2HeyI{t) 

R(t) = -H,{t)I{t) 

i{t) = He{t)R{t)-H,y(2P(t)-l) (3.6) 

where P = \ve\'^ is the Ve survival probability, R = Re(z/*i/e) and / = 
Im(z/*z/e). The corresponding initial conditions are: 

P(to) = l, /(to) = 0, R{to) = 0. (3.7) 

Defining 

R{t)±il{t) = e^'^*o^'=(*'^^*^Z±(t) 

Z±{t) = Ti f H,y{2P{t,)-l)e'^'-f^o''^^''^''''dh, (3.8) 

J to 

we can express the auxiliary functions R(t) and I(t) as: 

R{t) = f H,y{2P{U)-l)sin{f H,{t2)dt2)dtr (3.9) 

Jto Jti 

lit) = - f H,y{2P{ti)-l)cos{f He{t2)dt2)dU. (3.10) 

Jto Jti 

^The radial dependence of the solar matter density is understood as a time dependence 
since neutrinos are relativistic. 

•^In the Hamiltonian matrix, a term proportional to the identity has been removed. 



After substituting the Eqs. ( |3.9| ) and ( |3.10|) in the r.h.s of ( |3.6| ), we can 

average over the random density distribution, taking into account that for 
the random component we have: 



= 0, {A,y{t)A,y{h)) = 2K5{t - ti), (3.11) 
where the quantity n is defined as 

<t) = {Al{t))L, = \v^^{t){e)Lo- (3.12) 
At this point all we need are the following averaged products 

(^,(t)i?(t)) = -K(t)(/(t)), (ie,(t)/(t)) = <t){R{t)). (3.13) 



These are derived from Eqs. ( p.9| ) and ( |3.10| ) taking into account also 



and are correct up to 0{k). In terms of the averaged quantities defined 
as {Pit)) = P(U), {R{t)) = 7^(U), {I{t)) = J(U), we can write the noise- 
averaged variant of the set ( p.6|) as: 



Vit) = 2H,yl{U) 

n{t) = -2Aey{t)I{U) - G/^(U)7^(U) 

i{t) = 2A^y{t)n{u) - eK{u)l{u) -n^^{eV{u) - oo). (3.14) 

This system of equations explicitly exhibits the noise parameter k, 0. Elim- 
inating X and TZ from Eq. (|3.14| ) we can obtain the following third order 
differential equation for the averaged conversion probability V: 

Mt)Aey{t) - A,y{t) -^{t)+ U^,{t)A,y{t)+2A,y{t)k{t) 



dt 



dV / 



dt 

-2H%{A,y{t)-2K{t)A,y{t)), (3.15) 

where the frequency Uq familiar from the MSW effect is given as 

u;m=^{A%it) + Hl), (3.16) 

and the initial conditions become: 

P(U)=oo, PiU,) = f, #(U) = GHf^. (3.17) 

Let us notice that in the absence of noise {k = 0) the Eq. ( |3.15| ) reduces to 
the well known MSW equation (cfr. with (Eq. (2.23) of the first paper in 
Ref. lH) with the change H = 2H^y, H = 2A^y). 

In order to gain some more insight on the present picture let us note that 
the MSW resonance condition, i.e. Aey{t) = Vey{t) — Am"^ cos29/2E = 0, 
remains unchanged, due to the random nature of the matter perturbations. 
In other words, the fact that the noise is a second order effect (see eq. (|3 . 1 1|) ) 
means that it can only be seen in the conversion probability. In order to 
ensure that the correlation length Lq is smaller than the neutrino wave length 

^These equations are equivalent to those obtained in Ref. in terms of the variables 
x = 27^, y = -21 and r = 2V - oo. 



in the Sun, as required by the condition ( |2.7| ), we choose to adjust Lq as 
follows: 

Lo = 0.1 X (A™) = 0.1 X — . (3.18) 

In order to get a feeling for the importance of the noise term in the system 
( |3.14D , note that the noise parameter k in Eq. ( |3.12| ) is always smaller than 
Aey{t), for < few %, except at the resonance region. As a result, the density 
perturbation can have its maximal effect just at the resonance. However, this 
is not enough for the noise to give rise to sizeable effects. Since the noise 



term gives rise to a damping term in the system (|3.14|) , it follows that the 
corresponding noise length scale 1/k be much smaller than the thickness of 
the resonance layer Ar. In other words, it is also necessary that the following 
adiabaticity condition 

ttr = Ar Kres > 1, (3.19) 

is satisfied. This condition is analogous to the standard MSW adiabaticity 
condition > 1 where = Ar/ {Xm)res is the standard adiabaticity param- 



eter at resonance |T^. One can show that the two adiabaticity parameters 
are related as 

^2 Am^sm^2eRo 

ttr ~ «r n — ttr = 7^ — • (3.20) 

tan2 2^' 47rEcos2^ ^ ' 

For the range of parameters we are considering, ^ ~ 10~^ and tan^ 26 > 
10^^ — 10^^, and due to the restriction in the r.h.s of ( p.7|) one can estimate 
that Or < ar- Moreover, the relation a,,. < a,, can be rewritten as K^es < SHr, 
where SHr is the level splitting between the energies of the neutrino mass 
eigenstates at resonance. This shows that the energy perturbation induced 
by the matter fluctuations is not enough to cause the level crossing (even at 
the resonance) . In other words, it never violates the MSW adiabaticity 
condition f\ 

From Eq. (|3.20D it follows that, in the adiabatic regime > 1, the 



effect of the noise is larger the smaller the mixing angle value. Furthermore, 
as already noted above, Eq. ( p. 201) implies that the MSW non-adiabaticity 
ar < 1 is always transmitted to < 1- As a result, under our assumptions 
the fluctuations are expected to be ineffective in the non-adiabatic MSW 
regime. 



4. MSW Effect in Noisy Solar Matter 

In this section we study the impact that random perturbations in the so- 
lar matter density can have upon the MSW solution to the solar neutrino 
problem. We will consider both the active to active and active to sterile neu- 
trino conversion. For definiteness we will take as our reference SSM the most 
recent Bahcall-Pinsonneault (BP95) model with helium and heavy element 
diffusion, as given in the last paper in Ref . . From there we will take both 
the electron (neutron) density as well as the neutrino energy spectra and 

^This is opposite to the case of a local density jump as discussed by Krastev and 
Smirnov in the second paper in Ref. where larger values of 6p could break MSW 
adiabaticity. 



detection cross sections. Using these as input, we have solved numerically 
the coupled differential equations in ( ^.141 ) for the z/g survival probability 0. 

In order to get some preliminary insight on the effect of the density noise, 
in Fig.l we plot P as a function of E/ /S.m? for different values of the noise 
parameter ^. For comparison, the standard MSW case ^ = is also shown 
(lower solid curve). We take this the reference situation with which 

all others with non-vanishing ^ are compared. 

One concludes that in both cases of small and large mixing (Fig. la 
and Fig. lb, respectively), the effect of the matter density noise is to weaken 
the MSW suppression in the adiabatic regime (see dotted and dashed curves) 
with negligible effect in the non-adiabatic region, in complete agreement with 
the results of Ref . [|1^ . The relative increase of the survival probability V is 
larger for the case of small mixing (Fig. la) as already guessed on the basis 
of Eq. ( |3.2CI|) . One sees that the enhancement of the survival probability can 
easily reach 20% for ^ values as small as 4%. From these figures one can 
already infer that for the relevant Am? ~ 10"^ eV^ the intermediate energy 
neutrinos (like ^Be neutrinos) are the ones most likely to be affected by the 
matter noise. 

Note that the "white noise" -type density fluctuations we consider here 
cannot lead to any parametric enhancement of the survival probability 



of the type discussed in Ref. ^\ with a sinusoidal density perturbation. 



In contrast to that case, the effect of random perturbations is smooth, as 
suggested by the fact that the noise parameter k plays the role of a friction 
term in Eq. ( 3.14 ). 

Moreover, one can see from the figures that for the value E//S.m? ~ 
6.7 X 10^ eV~^ cos 29, required in order for the neutrinos to undergo resonant 
conversion just at the solar centre r = 0, the survival probability remains 
equal to 0.5 irrespective of the ^ values. The presence of this "fixed point" 
is easily understood: for such E / Aw? value Q the neutrino state ji/e >= 
v^d'^im > +|i^2m >)/2 is produced at its resonance point and ^ i^2m 
transitions between matter eigenstates occur at the same rate. This case of 
coincidence of neutrino production point with its resonance point is the only 
one for which the effect of the matter noise is strictly absent, even if the 
adiabaticity condition holds. 



5. Comparison with Solar Neutrino Experi- 
ments 

As seen above, there can be a substantial effect of matter noise fluctuations 
on the neutrino conversion probabilities. It is therefore important to analyse 
the possible impact of this scenario in the determination of solar neutrino 
parameters from the experimental data. The most recent averaged data of 

^For simplicity and CPU economy we have not included throughout our analysis the 
neutrino production distributions in the Sun. 

^Note that the matter mixing angle is understood as the mixing angle that diagonalises 
the Hamiltonian Eq. (|]^) at each time t. It is written as s\T?20m = '^H'^x/^o- ^t 
resonance ojq = 4iJg^, so that sin^ 29 „i — 1. 



the chlorine [0], galhum [Q,^ and Kamiokande experiments are: 



i?^7 = (2.55±0.25)SNU, = (74±8)SNU, R^!^ = (0.44±0.06)i?ff ^ 

(5.1) 

where R^a^^ is the BP95 SSM prediction. For the galhum result we have 
taken the of GALLEX i?Ga = (77±8±5)SNU § and SAGE i^^T = (69±11± 
6)SNU p measurements. The detection rates in the chlorine and gallium 
experiments are given as 

Rci,Ga = f dEa{E)V{£)Y.cl>^{8), (5.2) 

> 

where the sum is understood over the v source contributions {i =^Be, ^B...) 
and (t{E) are the corresponding neutrino cross sections. For the Kamiokande 
experiment, the detection rate is 

Rk= I dE\a,XE)V{£) + a,^{£){^-V{E))\(t>B{S), (5.3) 

JTh L ^ 

where a^^ (E) and 0"^^ (E) (x = /i, r) are the u^ — e and — e elastic scat- 
tering cross sections, respectively, and 'Th' stands for the detection energy 
threshold. In the case of sterile conversion a,y^ = 0. 

In Fig. 2 we show the iso-signal contours (within 2a standard deviations) 
for each experiment and for different ^ values. These plots demonstrate 
that the horizontal adiabatic lines are the ones mostly affected by the noise 
fluctuations. Indeed, the larger the ^ value, the greater the suppression of 
the neutrino conversion and, as a result, the larger the shift of this horizontal 
branch towards smaller Am^ values. These lower Am^ values allow neutrinos 
of lower energy to be involved in the adiabatic conversion (since the resonance 
matter density is proportional to Am^/E) so as to compensate for the effect 
of the matter noise. Notice also that, because of this downward shift in 
Am^, the minimum allowed values of sin^ 29 for each experiment becomes 
larger in order to preserve adiabaticity (see Eq. ( |3.2U| )). Also the diagonal 



(so-called "non-adiabatic") and the upper portion of the vertical (so-called 
large mixing) branches of the MSW plot are modified by the effect of noise. 
The diagonal lines are deformed mostly in the upper-left part, due to a shift 
of the kink towards larger values of the mixing angle. In contrast, they are 
less affected for sin^ 26' > 10~^. This follows from the fact that for larger 
mixing the noise adiabaticity is lost (see Eq. (|3.2CI| )). Thus we find that the 
expectations derived on the basis of our discussion in the previous section 
are confirmed. 

Comparing the allowed regions of all experiments shown in Fig. 2a, 2b 
and 2c for the or neutrino conversion case, one concludes that 

the overlapping area is not substantially changed in the small mixing branch, 
whereas it increases in the adiabatic one, for large mixing sin^26' > 0.3 and 
lO-^eV^ < Am2 < IQ-W. 

We now turn to the case of '^s sterile resonant transitions. Here the 
z/e survival probability is not substantially changed with respect to the active 
neutrino conversion case, since the solar neutron contribution in the matter 
potential is rather small, compared to that of the electrons. As a result the 
signal expected in radiochemical experiments is rather insensitive to whether 
the converted neutrino is active or sterile. Thus we focus on the Kamiokande 



experiment. In Fig. 2d we show the iso-signal contours for the case of sterile 
neutrino conversion. One can see from the figure that in the sterile case, 
irrespective of the assumed level of noise, the vertical large mixing branch 
gets thinner and closer to the maximal mixing region sin^2^ = 1. This is 
required, of course, in order to increase the contribution to the signal which 
is now lost when compared to the active One can also see that the 

noise has the same qualitative features as in the case of active conversions, 
mostly affecting the horizontal adiabatic region of larger Am^ . 

In order to determine the solar neutrino parameters Am^ and sin^ 29 we 
now proceed to perform a analysis for the present experimental data. For 
simplicity we neglect for the moment the theoretical uncertainties. 

The results of our fit are shown in Fig. 3, where the 90% confidence 
level (C.L) areas are drawn for different ^ values (see also Table 1). From 
Fig. 3a one can observe the modifications in the small mixing region caused 
by the noise in the case of active neutrino conversion. One sees that there 
is a slight shift of Am^ towards lower values and a larger shift of sin^ 29 
towards larger values. For example the allowed region for the mixing angle 
covers the range 4 x 10^"^ < sin^ 2^^ < 8 x 10^"^ obtained for = becomes 
8 X 10"^ < sin^26' < 2 x 10"^ for ^ = 8%. The corresponding allowed 
Am^ range is 2.5 x 10~^ < Am^ < 9 x 10~^ eV^ to be compared with 
5 X 10~^ < Am^ < 1.2 X 10~^ eV^ in the noiseless case. The large mixing 
area is less stable, exhibiting a tendency to shift towards smaller Am^ and 
sin^2^^. For example, if we take ^ = 8%, for the sake of argument, we find 
that the small mixing region is much more stable than the large mixing one, 
even for such a relatively large value of the noise. 

As for the value of the minimal x^, the presence of the matter density 
noise makes the data fit a little poorer: Xmin = 0.1 for ^ = 0, Xmin — 0-8 
for ^ = 4% and Xmin — 2-1 for ,^ = 8%. Also the best fit points where 
the Xmin achieved change slightly: the value of Am^ ~ 6 ^ 7 x 10~^eV^ 
is almost unchanged, while the value of the mixing angle gets larger with 
respect to the noiseless case. For example, sin^ 2^^ = 6 x 10~^ for ^ = while 
sin^ 2^ = 8 X 10-3 for ^ = 4% and sin^ 29 = IQ-^ for ^ = 8%. The strong 
''Be neutrino suppression, characteristic of the MSW effect, is reduced by 
the presence of matter noise (see Fig. 1). As a result, the confiict between 
chlorine and Kamiokande data is exacerbated and the data fit gets worse. In 
any case our results for Xmin (see Table 1) indicate that the MSW scenario 
still provides a good fit of the totality of solar neutrino data, even in the 
presence of matter fiuctuations, as long as ^ < 8% or so. 

As for the large mixing solution, although the Xmin value is not substan- 
tially changed with respect to the noiseless case (see Table lA), we find that 
it acquires an increased statistical significance with respect to the correspond- 
ing region of the noiseless case. Our results show that the large mixing solu- 
tion gets wider than in the noiseless case. For example the smallest allowed 
sin^ 29 value shifts from sin^ 29 ~ 0.4 for C = 4% down to sin^ 29 ^ 7 x 10"^ 
for C = 8%. 

As for the best fit points we find that Am^ ~ 10~^eV^ is almost un- 
changed, whereas the best value of the mixing angle decreases from sin^ 29 = 
0.67 for the noiseless case down to sin^2^ = 0.27 for ^ = 8%. Note that 
the possibility of lowering the mixing angle value sin^ 29 characterising the 
large mixing MSW solution in the presence of noise may eliminate the super- 



nova argument given in Ref. p5i against such solution. In agreement with 



Ref. |T3[, we find that in our fit this region appears aheady at the 80% CL. 
in the ^ = case f\. 

We now turn to the case of sterile solar neutrino conversions. We find 
that the data fit is worse {xLin — 1) than for the active case (see Fig. 3b 
and Table IB) and it excludes, even at 95% C.L., the large mixing region 



(in the noiseless case). This is in agreement with previous analyses |15|,p!3 
However, the presence of matter density noise may restore this region. For 
example for ^ = 8%, although the data fit is much worse than in the ^ = 
case, the large mixing region appears at the 90% CL. We may note in this 
context that the indicated range for the mixing angle is not in confiict with 
the primordial helium abundance constraints p6[| . 

So far in our analysis we have neglected SSM theoretical uncertainties, 
and worked entirely within the BP95 model p. One way to account for 
these uncertainties would be to allow the solar neutrino fiuxes to vary as 
suggested in Ref. |]16|,|12|. However one can get an idea (even if partial) 
of these uncertainties by simply repeating the data fit assuming the SSM 
of Turck-Chieze et al. (TCL) 0. For our purposes the main difference 
between this model and the BP95 model is that it predicts a lower flux. 
The comparison of the allowed parameter regions obtained in the framework 
of the TCL model. Fig. 3 (c,d), with those obtained using the BP95 model, 
Fig. 3 (a,b), shows that the general features of the effect of the noise are 
maintained. In particular, our results once again establish the fact that the 
indicated Am? range for the small mixing MSW solution is fairly stable, as 
long as the assumed noise level is not too large. Note also from the flgures 
that, even though the effect of the noise is to lower the Am? range for the 
large mixing solution, the region obtained (e.g. for ^ = 8%) lies higher than 
the corresponding range for the BP95 model. 



6. Implications for Future Experiments 

Up to now we have discussed the possible consequences of the presence of 
matter fluctuations for the ongoing solar neutrino experiments. We now turn 
our attention to the possibility of probing the level of matter noise in the Sun 
in the next generation of solar neutrino experiments. 

As we have seen the ^Be neutrinos are the component of the solar neutrino 
spectrum which is most affected by the presence of matter noise. Therefore 
the future Borexino experiment, aimed to detect the ^Be neutrino flux p7 



through the elastic v — e scattering should be an ideal tool for studying the 
solar matter fluctuations. 



In Ref. IjT^ it was shown that in the relevant (noiseless) MSW parameter 
region the Borexino signal cannot be sharply predicted. This is illustrated 
in Fig. 4a, where we display the Borexino signal in the Am^ - sin^ 29 plane, 
expressed in units of the expected SSM rate, i.e. Zse = R^bI'^ / RbV^ ■ one 
can see, the allowed range for the signal in this case lies anywhere between 
0.2 to 0.7 of the SSM prediction. In Fig. 4b, we show the corresponding 
beryllium line predictions for the case of noisy MSW, assuming ^ = 4%. 



^This result may be underestimated since the earth regeneration effect has not been 
included. 



We see that the presence of matter noise strongly modifies the picture: the 
minimal allowed value for Zse now becomes higher, Zse > 0.37. Therefore 
if the Borexino experiment detects a small signal, Zse ^ 0.3 (with sufficient 
accuracy) this will imply that a 4% level of matter fiuctuations in the central 
region of the Sun is rather unlikely to be present if the MSW mechanism is 
responsible for the explanation of the solar neutrino deficit |. 

Note, on the other hand, that if a higher value Z^e ^ 0.5 would be found 
experimentally, this would be incompatible with the small mixing MSW so- 
lution with noise at the ^ = 4% level. However, this higher signal could 
be consistent with the both the large mixing MSW solution as well as the 
noiseless small angle MSW solution. On the other hand, if the noise level is 
higher, ^ = 8%, the allowed Zse range narrows down to values between 0.5 
to 0.65. 

Let us turn to the case sterile resonant conversion in the noisy MSW 
effect. Let us imagine that future large detectors such as Super-Kamiokande 
and/or the Sudbury Neutrino Observatory (SNO) establish through, e.g. 
the measurement of the charged to neutral current ratio, that the deficit of 
solar neutrinos is due to the v^. — > Vg resonant conversion. In this case, the 
minimum signal expected in Borexino is very small Zse ~ 0.02 for ^ = (see 
Fig. 4c). On the other hand in the noisy case with ^ = 4%, the minimum 
expected Borexino signal is 10 times higher than in the noiseless case, so that 
if Borexino detects a rate Zse ^0.1 (see Fig. 4d) this would again exclude 
noise levels above 4%. 

7. Discussion and Conclusions 

We have presented a comprehensive study of the effects of the matter den- 
sity noise upon the MSW solution to the solar neutrino problem. We have 
adopted the wave function Schrodinger formalism to re-write the correspond- 
ing MSW evolution equations for the neutrino survival probabilities. The 
fluctuations weaken the efficiency of the MSW suppression in the adiabatic 
regime, whereas they are much less effective in the non-adiabatic regime. In 
our data flt we have shown that the MSW solution still exists for realistic 
levels of matter density noise ^ < 8%. However, our analysis has shown 
that the quality of the flt gets a little worse if these noisy matter pertur- 
bations are present. In any case the mass range determined from our flt 
for the small mixing MSW solution 4 x 10~^eV^ < Am^ < 10~^eV^ is rela- 
tively stable at 90% C.L., whereas the mixing angle determination appears 
more sensitive to the assumed level of fluctuations, and shifts sin^ 2Q towards 
larger values up to 10~^. These trends also hold for the case of sterile solar 
neutrino conversion. In the latter case we have found that in the presence 
of solar density noise the large mixing region gets somewhat improved sta- 
tistical signiflcance when compared with the noiseless case. However, it is 
remains highly disfavoured with respect to the small mixing MSW solution. 

We have also explored the potential of the Borexino experiment to "test" 
the level of matter density fluctuations in the solar interior through the mea- 

^In principle any value of Zbc is also compatible with the just-so oscillation scenario 
p^p^ , but here the strong seasonal ''Be and pep signal variations, would help to distin- 
guish from the MSW case. 



surement of the ^Be neutrino flux, as depicted in Fig. 4. 

Finally, we note that in our analysis we have neglected the details of 
the neutrino production distribution as a function of the distance to the 
solar centre. It is well known that this affects mainly the low energy pp 
neutrinos [|1^]. As a result, the iso-signal curves we have obtained for the 
gallium experiments are somewhat less reliable in the position of the kink 
corresponding to Am^ ^ 2 x 10~^eV^ marking the onset of pp neutrino 
suppression and lying on the gallium non-adiabatic branch. However, this 
does not substantially affect the determination of the relevant regions where 
all solar neutrino data are explained through the MSW effect. This includes 
both small as well as large mixing MSW solutions. 
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A: active 


C = o 


C = 2% 


C = 4% 


C = 8% 


small 9 












0.10 


0.23 


0.80 


2.1 


Am2(10-^eV2) 


0.68 


0.71 


0.65 


0.61 


sin^ 26 


6.2 X 10-3 


7.3 X 10-3 


7.5 X 10-3 


10-2 



large 6 












3.3 


2.9 


3.0 


3.2 


Am2(10-5eV2) 


2.7 


2.4 


2.0 


1.2 


sin^ 26 


0.67 


0.69 


0.57 


0.27 



B: sterile 


( = 


C = 2% 


e = 4% 


C = 8% 


small 6 












1.0 


1.9 


3.6 


8.9 


Am2(10-^eV2) 


0.53 


0.50 


0.49 


0.40 


sin2 29 


7.5 X 10-3 


7.5 X 10-3 


9.0 X 10-3 


1.3 X 10-2 



large 9 










X^ ■ 


10 


11 


11 


11 


Am2(10-^eV2) 


1.4 


1.2 


1.6 


1.0 


sin2 29 


0.83 


0.83 


0.69 


0.39 



Table 1. The values of Xmin^ ^'^^ ^ degree of freedom, and the corresponding 

best fit Am2 and sin2 29 parameters in the small and large mixing regions, for 
different values of ^. Tables A and B arc for the active-active and active-sterile 
conversion respectively, using the latest 1995 Bahcall-Pinsonneault (BP95) model. 



A: active 
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e = 2% 


C = 4% 


C = 8% 
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0.10 


0.46 


1.1 


3.0 


Am2(10-5eV2) 


0.64 


0.66 


0.62 


0.66 


sin^ 261 
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3.6 X 10-3 
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5.2 
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4.3 


Am2(10-5eV2) 


17 


16 
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3.8 


sin^ 29 
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0.83 
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0.69 



B: sterile 
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C = 8% 


small 9 












0.69 


2.1 


2.3 
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Am2(10-^eV2) 


0.51 


0.45 


0.48 


0.57 


sin^ 29 


4.3 X 10-3 


4.2 X 10-3 


5.2 X 10-3 


6.3 X 10-3 
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9.3 


9.3 


9.2 


8.3 


Am2(10-^eV2) 


16 


15 


14 


3.3 


sin^ 29 


0.83 


0.83 


0.83 


0.69 



Table 2. The same as for Table 1, but using the Turck-Chieze et al. (TCL) 
SSM. 



Figure Captions 



Fig. 1. 

Averaged solar neutrino survival probability V versus Ej [^vr? for small mixing (a: 
sin^ 10 = 10~^) and large mixing (b: sin^ 29 = 0.7). The solid, dotted, dashed and 
dot-dashed curves correspond to noise levels ^ = 0, 2%, 4% and 8%, respectively. 

Fig. 2. 

Iso-rate contours for the chlorine (a) , gallium (b) and v — e scattering (c) experi- 
ments for the case of active neutrino conversion, ^h,t- The threshold energy 
for the recoil electron detection is 7.5 MeV. For the radiochemical experiments the 
results are in SNU, whereas for the v — e scattering experiment these are given in 
units of the BP95 SSM prediction. The contours delimit the 2a allowed regions. 
The solid, dashed and dotted curves correspond to = 0, 4% and 8%, respectively. 
Fig. 2d gives the same iso-rate contours for the case of sterile conversion 
for the u — e scattering experiment. 

Fig. 3. 

The 90% C.L. allowed regions (given by the condition < Xmm + 4.61) for the 
active (a and c) conversion and for the sterile (b and d) conversion. For Fig. 
3a and 3b, the C.L. allowed regions are obtained using the most recent Bahcall 
and Pinsonncault (BP95) SSM; for Fig. 3c and 3d we used the Turck-Chieze and 
Lopes (TCL) SSM. In Fig. 3a, 3b, 3c and 3d the solid, dot, dash, dot-dash curves 
correspond to the cases ^ = 0,2%, 4% and 8%, respectively. The Xmin ^^"^ 
corresponding (Am^, sin^ 29) best fit points are given in Table 1 and Table 2. 

Fig. 4. 

The iso-signal contours of the ratio Zse = R^Be'^ / ■^Be'^ (figures at the curves) 
in the u — e scattering Borexino detector (solid lines) . The threshold energy for 
the recoil electron detection is 0.25 MeV. The 90% C.L. allowed regions (dotted 
lines) and the corresponding best fit points (diamonds) are also superimposed, as 
determined by the present experimental data and using the BP95 SSM. 

The case of active resonant conversion is presented in Fig. 4a and Fig. 4b 
for ^ = and ^ = 4%, respectively. Analogously, Fig. 4c (^ = 0) and Fig. 4d 
(C = 4%) refer to the case of sterile neutrino resonant conversion. 
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